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ABSTRACT 

Aims. We aim at measuring mass-loss rates and the luminosities of a statistically large sample of Galactic bulge stars at several 
galactocentric radii. The sensitivity of previous infrared surveys of the bulge has been rather limited, thus fundamental questions 
for late stellar evolution, such as the stage at which substantial mass-loss begins on the red giant branch and its dependence on 
fundamental stellar properties, remain unanswered. We aim at providing evidence and answers to these questions. 
Methods. To this end, we observed seven 15x15 arcmin 2 fields in the nuclear bulge and its vicinity with unprecedented sensitivity 
using the IRAC and MIPS imaging instruments on-board the Spitzer Space Telescope. In each of the fields, tens of thousands of point 
sources were detected. 

Results. In the first paper based on this data set, we present the observations, data reduction, the final catalogue of sources, and a 
detailed comparison to previous mid-IR surveys of the Galactic bulge, as well as to theoretical isochrones. We find in general good 
agreement with other surveys and the isochrones, supporting the high quality of our catalogue. 

Key words. Galaxy: bulge — Galaxy: stellar content — Infrared: stars — Stars: late-type — Stars: mass-loss — Stars: AGB and 
post-AGB 



1. Introduction 

The Galactic bulge (GB), an important dynamical and morpho- 
logical component of our Galaxy, offers an environment distinct 
from the Galactic disk for study of stellar populations, stellar 
evolution, and the mass-loss processes that accompany and, in 
the end, control the last. Understanding and calibrating the phys- 
ical processes whereby mass ejected by evolved stars into the 
bulge environment is recycled back into new generations of stars 
requires a statistical knowledge of mass loss as a function of fun- 
damental stellar parameters in this region. Because of the limited 
sensitivity of previous surveys of the bulge, fundamental ques- 
tions for late stellar evolution, such as the stage at which substan- 
tial mass-loss begins on the red giant branch (RGB), and its de- 
pendence on fundamental stellar properties, remain unanswered. 
The GB is an ideal laboratory for addressing these issues, provid- 
ing a very large sample of stars at an almost identical distance. 

We therefore observed seven 15 x 15 arcmin 2 fields that sam- 
ple a range of distances from the Galactic centre with unprece- 
dented sensitivity using the Infrared Array Camera (IRAC; Fazio 
et al.|2004] l and the Multiband Imaging Photometer for Spitzer 



(MIPS; Ri eke et al.|2004| l, the imaging instruments on-board the 
Spitzer Space Telescope (Werner et al. 2004 1, in order to de- 



* Table 5 as well as the fits mosaics are only available in elec- 
tronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-t-A/ 



termine mass-loss rates and luminosities of a statistically large 
sample of stars at several galactocentric radii. These data en- 
able us to detect stars with very low mass-loss rates through 
their infrared excess, determine the dependence of the mass- 
loss rate on luminosity and effective temperature along the giant 
branches, and conduct a census of mass-losing stars at different 
rates. The observations, together with existing studies that probe 
higher mass-loss rate stars, will enable us to infer the total rate 
of mass loss in the bulge, a key input to evolutionary models of 
the bulge. The data have already led to the discovery of mid-IR 
logP vs. magnitude relations ( Glass et al.||2009| l. In this paper 
we present the observations, the data reduction, and the source 
catalogue. We also compare the data to previous mid-IR surveys 
of the bulge and to theoretical isochrones. 

The outline of this paper is as follows. In ^2] we describe 
the observations, the data reduction, the different steps of point 
source extraction, and how the catalogues were created. We then 
proceed in ^3] with checking our photometric data against other 
mid-IR surveys of the Galactic bulge. Further checks are pre- 
sented in ^4] where colour-magnitude diagrams (CMDs) are 
compared to theoretical isochrones. Finally, ^5] draws conclu- 
sions on the catalogues and the data quality. 
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2. Observations and Data Reduction 



2.1. Field selection and observations 

The locations of the observed fields were chosen to sample the 
bulge on a variety of scales, to measure how the mix of stellar 
populations varies with Galactic latitude. They were also chosen 
to avoid the relative intense, saturating emission from near the 
Galactic plane. Our innermost fields are within the central stel- 
lar cusp, which presumably contains stars of various ages (c.f. 
Blum et al. 2003 ; Figer et al. 2004). Because the stars within this 
domain are believed to have formed within the central molecu- 
lar zone and then to have diffused into an increasingly thicker 
distribution as a result of scattering off molecular clouds ( |Kim| 
& Morris 2001 ), we expect a vertical segregation of stellar ages. 
Thus, inasmuch as the luminosities and mass-loss rates of red 
giants and asymptotic giant branch (AGB) stars depend on their 
masses, hence their ages, the radial distributions of the different 
kinds of evolved stellar objects can be used to model the star 
formation and dynamical history of this region. 

Two fields were selected to sample the nuclear bulge at 
(l,b) = (0?00,-l°00) and (l,b) = (0?63,-0?36), namely Nl 
and N2, both of which were observed by ISOGAL (project 
for imaging part of the Galaxy using ISO, the Infrared Space 
Observatory; Omont et al. 2003) at 7 and 15 p.m. Five fields 
were selected beyond the nuclear bulge in areas where they over- 
lap the Optical Gravitational Lensing Experiment III (OGLE-III; 
Udalski et al.|2000) micro-lensing survey. Four of these are lo- 



cated below the Galactic plane along a radial vector that sub- 
tends the minor axis of the Galaxy at an angle of about 13°5, 
and terminates at the well-studied field in Baade's window at 
(I, b) = (V.03, -3?83), c ontaining the globular cluster NGC 6522 
(e.g. Glass et aI7||1999) >. The three fields inside this are located 
at (l,b) = ((f.30,-rA2), (0!56, -T.23), and (0?76, -3?07). The 
fifth field is positioned above the plane at (l,b) = (2?87,0?35). 
All fields are approximately rectangular in right ascension (RA) 
and declination (Dec) coordinates. Figure [T] shows the location 
of the seven observed fields with respect to the Galactic cen- 
tre, and Table[T] summarises some of their main characteristics. 
The range in RA and Dec (J2000) given in Table [T] refers to 
where there is full overlap between MIPS and all four IRAC 
bands. Particularly in IRAC we do have some coverage outside 
the given range. 

Observations of the GB fields were performed using the 
IRAC instrument in all four channels at 3.6, 4.5, 5.8, and 8.0 pm, 
and the MIPS instrument in the 24 pm channel on board the 
Spitzer Space Telescope within the programme 2345. Bulge field 
N 2 was observed only with MIPS in our programme, because 
this field was covered by IRAC observations within the gen- 
eral observer (GO) programme 3677 (principal investigator: S. 
Stolovy; Ramf rezet al.|2 008). The IRAC observations were car- 
ried out on 2005 March 30 and 31, and the MIPS observations 
between 2005 April 08 and 13. 

With IRAC, observations were done in the full-array read 
out-mode, one frame per pointing with 2 s integration time per 
frame - the shortest full-array integration time was used to min- 
imise the effects of saturation and resulting latency problems. 
The mapping was done in a 6 x 6 rectangular grid with step size 
260", with five dither positions, and a medium scale factor, giv- 
ing a total exposure time of 10 s per pixel. The MIPS observa- 
tions were obtained using the photometry/raster mode with 10 s 
integration time and full-array read-out mode, 3x3 rectangular 
grid and two cycles, giving a total integration time of 331 s per 
pixel. 
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Fig. 1. Map of the seven fields observed in this campaign with 
the labels used throughout this paper. 



The fields of view of IRAC channels 1 and 3, 2 and 4, as well 
as MIPS did not fully overlap. Nevertheless, the MIPS field of 
view is fully contained in all four IRAC fields of view. 

2.2. Data reduction 

We corrected the IRAC basic calibrated data (bed) files to mit- 
igate artifacts such as muxbleed or column pulldown with the 
tools provided by S. Carey at the Spitzer Science Center (SSC). 
Post-bed processing was then conducted on the corrected bed 
files using the MOsaicker and Point source Extractor (MOPEX) 
software and its subsystem APEX (version 18.2.0, Markovoz & 
|Marleau|2005) . 



2.2.1 . Step 1 : Point source detection and extraction 

Before further processing, mosaics were created by MOPEX 
for each field and channel from the corrected bed frames. In 
this step, the pipeline interpolates the input images onto the 
output grid, taking geometric distortion into account. An out- 
lier detection scheme flags bad pixels and any pixels affected 
by cosmic ray hits or moving objects, and these pixels are re- 
computed before co-addition. Finally, the interpolated images 
are co-added to one mosaic image. After mosaicking, APEX de- 
termines the background by calculating the median in a 45 X 45 
pixel box around each pixel, and subtracts it from the image. 
These background-subtracted images were used in the detec- 
tion step. Then, background fluctuations in the images were esti- 
mated and noise images derived, which were used for signal-to- 
noise ratio (SNR) estimation and the generation of point source 
probability images. With those, the detection table was com- 
piled. We chose a detection threshold of 2.3cr above the back- 
ground. We used both the point response function (PRF) fitting 
capability and the aperture photometry functionality with circu- 
lar apertures with radii of 2, 3, and 4 pixels. Larger apertures 
were not applicable owing to the crowding in the fields. The tiles 
used for PRF fitting also had to be chosen to be very small (3x3 
pixels). We used the most up-to-date mean PRFs as provided by 
the SSC, and set the PRF normalisation radius accordingly. 

The measured x 1 values in our data seem to be very high, 
even for successfully fitted sources, e.g. in the least-crowded 
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Table 1. Main characteristics of the seven fields towards the GB observed with Spitzer as part of programme 2345. 
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Field name 


RA range (h m s) 


Dec range (° ' ") 


1 centre 


b centre 


Date IRAC obs 


Date MIPS obs 


Bulge 1 


17 51 14.2 . 


.17 52 36.7 


-29 15 22 . 


.-29 34 15 


0.30 


-1.42 


2005 03 30 


2005 04 09 


Bulge 2 


17 55 02.1 . 


. 17 56 24.2 


-29 26 25 . 


.-29 45 13 


0.56 


-2.23 


2005 03 31 


2005 04 10 


Bulge 3 


17 58 50.7. 


.18 00 13.3 


-29 41 08 . 


.-29 59 59 


0.76 


-3.07 


2005 03 30 


2005 04 10 


Bulge 4 


17 50 16.1 . 


. 17 51 36.4 


-26 08 56 . 


.-26 27 51 


2.87 


0.35 


2005 03 30 


2005 04 09 


Bulge N 1 


17 48 51.9. 


. 17 50 14.8 


-29 17 42 . 


. -29 36 42 


0.00 


-1.00 


2005 03 30 


2005 04 08 


Bulge N 2 


17 47 55.8. 


.17 49 16.3 


-28 28 41 . 


.-28 46 19 


0.63 


-0.36 




2005 04 13 


NGC 6522 


18 02 30.4 . 


.18 03 53.5 


-29 49 26 . 


. -30 08 24 


1.03 


-3.83 


2005 03 30 


2005 04 10 



field NGC 6522, they are between 2.5 and 3.3 in the IRAC chan- 
nels 1 and 2, between 1 .4 and 1 .8 in IRAC 3, and only in IRAC 4 
at the perfect value of 1 . This may stem from confusion noise, 
which is not included in the provided uncertainty images and/or 
from the high source density. Confusion noise is created by 
the amplitude variations from PRFs of closely spaced sources 
( |Hacking & Houck|1987"l|Rieke et al.|1995| ). Modelling the con- 
fusion noise has not improved the resulting x 2 values. This has 
consequences for the reliability of automatic de-blending proce- 
dures, which were not working properly in this field. We thus 
had to disable the active and passive de-blending capabilities of 
APEX (see Appendix [A). 

In all MIPS fields, we first detected bright sources, removed 
the Airy rings of these sources to do the source detection on the 
residual images, and performed the point source extraction on 
the original images. The southern half of the field N 2 has strong 
diffuse emission. In that area, the detection threshold is certainly 
higher than in areas with low diffuse background emission, as is 
also evident from the histograms of magnitudes in Appendix [B| 



2.2.2. Step 2: Assigning quality grades to the sources 

Three different effects can be seen when plotting the flux ob- 
tained by PRF fitting versus that found with aperture photome- 
try for all the extracted sources (Fig.|2|. The first effect is present 
when we use a large aperture with a radius of 3 or 4 pixels. At 
very low PRF fluxes, there are several sources that show dispro- 
portionately large aperture fluxes, i.e. are above the dashed line. 
These sources are faint sources close to a bright source. In case 
of the larger aperture, we therefore already get a significant con- 
tribution in the aperture from the bright source. The second effect 
is relevant for some bright sources, whose PRF fluxes are smaller 
than the aperture fluxes. These sources are super- saturated; i.e., 
the radial intensity profile has a dip in the centre and so it shows 
a double peak structure. In a few cases, depending on the depth 
of the central valley, these two peaks are identified as two indi- 
vidual sources that are each fitted with a fainter PRF; i.e., both 
sources lie above the dashed line. In other cases, the fitted PRF 
ends up with a central flux between that of the central valley and 
that of the wings, meaning PRF fitting underestimates the total 
flux, so these sources are also shifted to the left in Fig. [2] 

The SSC homepage provides values of the maximum flux 
of unsaturated point sources as a function of integration time. 
Table |2]presents the flux values interpolated to a frame exposure 
time of 1.2 s, along with a lower flux limit. This lower flux limit 
was adopted on subjective grounds for the quality grade labelling 
(see below). We thus have to be cautious with sources brighter 
than about 320 mJy in IRAC 1 . Consistent with this warning is 
a third slight effect that seems to set in at fluxes close to this 
value and affects all brighter sources: the PRF fluxes increase 
disproportionately compared to the aperture fluxes, leading to a 
bending to the right, away from the dashed line in Fig. [2] (top). 



Table 2. Flux limits applied for assigning quality grades. 



Band 


Lower limit (mJy) 


Saturation limit (mJy) 


IRAC 3.6/an 


40 


320 


IRAC 4.5 yum 


30 


330 


IRAC5.8yum 


20 


2300 


IRAC 8.0//m 


20 


1200 


MIPS 24 nm 


40 


220 



Aperture photometry thus measures the saturated plateau of the 
PRF, while PRF fitting seems to ignore the missing flux that 
would be present beyond the plateau and fits the wings. The indi- 
vidual source then lies below the dashed line, with a larger PRF 
flux than the aperture flux. 

If dividing the aperture flux by the linear least squares fit to 
sources with x 1 < 5, the fluxes measured with both ways are 
almost identical for many of the sources with^ 2 < 5. The slope 
of this fit is not the aperture correction factor, to be derived in 



Sect. 2.2.3 Except for the sources affected by the three afore- 
mentioned effects, the relative difference calculated as 



rel. difference 



f f 

r PRF _ I- 



corrected 
aperture 



(1) 



scatters around 7% (Fig. [2] bottom). The superscript "corrected" 
in Eq.[T] indicates that the aperture flux has been divided by the 
linear fit to sources with x 1 < 5. Sources that are within this 
range are considered as quality grade A. This test is performed 
first, therefore grade A sources span the whole range of fluxes, 
from faint to almost saturated sources. 

Faint sources with fluxes smaller than a certain lower limit 
(Table [2} and outside the 7% range are classified as quality 
grade B. Sources whose fluxes are above the saturation limit are 
classified as quality grade C. They would have to be treated dif- 
ferently with dedicated tools. However, since these sources are 
probably mostly foreground stars, they are not of prime interest 
to the science goals, thus we decided not to invest more time 
in their flux determination. The interested reader is invited to 
download and analyse the mosaics, which will be made avail- 
able as on-line material at Centre de Donnees astronomiques de 
Strasbourg (CDS). 

Sources that failed all of the above criteria are classified as 
quality grade D. These are blended sources for which the PRFs 
overlap such that even the aperture with the smallest radius in- 
cludes both peaks; i.e., no reliable flux measurement is available 
for these sources. This effect is the same as the first one men- 
tioned above related to the larger aperture with a radius of 3 pix- 
els or larger. 

The statistics of grades for all seven fields and all five bands 
are given in Table [3] The percentage given for grades A and B 
is relative to the total number of sources in the catalogue for the 
respective field that is given in the last column. The seventh col- 
umn gives the number of sources not detected in the respective 
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Table 3. Statistics of sources with their grades for all seven fields and five bands. 



Field 


Band 


Grade A 


Grade B 


Grade C 


Grade D 


Not detected 


total 


Bulge 1 


3.6 pm 


12446(18%) 


39177 (57%) 


183 


626 


16661 


69093 




4.5 /im 


16173 (23%) 


35598 (52%) 


85 


216 


17021 






5.8pm 


13839 (20%) 


21757 (31%) 


1 


169 


33327 






8.0pm 


12669(18%) 


17044 (25%) 


7 


96 


39277 






24 pm 


577 ( 1%) 


1197(2%) 


1 


2 


67316 




Bulge 2 


3.6pm 


14889 (25%) 


30525 (51%) 


197 


748 


13963 


60322 




4.5 pm 


17495 (29%) 


28190(47%) 


72 


492 


14073 






5.8pm 


13712 (23%) 


17900 (30%) 





214 


28496 






8.0 pm 


11960 (20%) 


15752 (26%) 


1 


152 


32457 






24 pm 


503 ( 1%) 


1119(2%) 


4 


4 


58692 




Bulge 3 


3.6pm 


13476 (26%) 


27927 (53%) 


131 


101 


11031 


52666 




4.5 pm 


16329 (31%) 


25159(48%) 


62 


64 


11052 






5.8pm 


12085 (23%) 


15880 (30%) 





35 


24666 






8.0 /an 


10712 (20%) 


13220 (25%) 


4 


24 


28706 






24 pm 


383 ( 1%) 


976 ( 2%) 


1 


2 


51304 




Bulge 4 


3.6pm 


12938 (22%) 


33046 (57%) 


136 


197 


12091 


58408 




4.5 pm 


15689 (27%) 


29961 (51%) 


75 


95 


12588 






5.8pm 


11311 (19%) 


18283 (31%) 


3 


177 


28634 






8.0 pm 


9239 (16%) 


13588 (23%) 


19 


98 


35464 






24 pm 


267 ( 0%) 


817 ( 1%) 


3 


4 


57317 




Bulge N 1 


3.6pm 


12553 (19%) 


35826 (55%) 


198 


724 


15534 


64835 




4.5 pm 


15402 (24%) 


34527 (53%) 


85 


296 


14525 






5.8pm 


11819(18%) 


21845 (34%) 


5 


520 


30646 






8.0//m 


12035 (19%) 


16565 (26%) 


11 


144 


36080 






24 pm 


467 ( 1%) 


1126(2%) 


3 


23 


63216 




Bulge N 2 


24 pm 


21 ( 5%) 


316(74%) 


17 


74 





428 


NGC 6522 


3.6pm 


12430 (26%) 


24924 (53%) 


79 


55 


9796 


47284 




4.5 pm 


15637 (33%) 


21263 (45%) 


32 


21 


10331 






5.8pm 


12011 (25%) 


11907 (25%) 





30 


23336 






8.0pm 


10263 (22%) 


10432 (22%) 





12 


26577 






24 pm 


431 ( 1%) 


1152(2%) 








45701 





channel, but detected in at least one of the other bands plus in Table 4. Aperture correction factors 



the reference catalogues (see Sect. 2.2.4 1. 



2.2.3. Step 3: Determination of the aperture corrections 

Since we do not use a larger aperture to estimate the background, 
which would give unreliable photometry in the crowded fields, 
but do aperture photometry on the background-subtracted im- 
ages instead, the aperture corrections are expected to be different 
from the standard values given in the instrument handbooks. 

To derive aperture corrections, we created point spread 
functions with stinytim2.0 (Krist 2006} for all five bands. We 
then inserted artificial sources with known flux in our mo- 
saicked, background-subtracted images and compared the result- 
ing fluxes given by APEX with the initial fluxes with which we 
had created the sources. The slope of a linear least squares fit to 
the initial vs. measured flux gives the aperture correction factors. 
The derived aperture corrections for the aperture with two pixels 
radius are given in Table |4] along with their statistical error as 
derived from the linear fit. Also given in the table are the values 
as recommended by the SSC, for a 2 pixel radius on source and 
a 10 - 20 pixel radius background annulus in the case of IRAC, 
and a 3'.'5 on source aperture without any background annulus 
in the case of MIPS. The SSC recommended values are not di- 
rectly comparable to our aperture correction factors, so they are 
only given for illustrative purposes. Despite the high value of 
the aperture correction factor for MIPS 24, we adopted the small 
aperture with 2 pixels radius because i) our observations are not 
photon noise limited, thus a larger aperture would not reduce 
the error estimate; ii) there is a "population" of blended sources 



Band 


Aperture correction 


SSC recommended 


IRAC 3.6pm 


1.2129±0.0008 


1.205 


IRAC 4.5 pm 


1.2730±0.0005 


1.221 


IRAC 5.8pm 


1.4693±0.0007 


1.363 


IRAC 8.0pm 


1.6556±0.0006 


1.571 


MIPS 24 pm 


1.5742±0.0007 


2.560 



whose flux increases disproportionally with growing aperture 
because more of the neighbouring source's flux is measured with 
larger aperture; iii) the SSC lists aperture corrections for an even 
smaller aperture, so our value is not an extreme. We also cor- 
rected the PRF fluxes with the photometric correction factors as 
recommended by the SSC. 

Besides the aperture correction and the photometric correc- 
tion for the PPvF fluxes, we did not apply any other correction 
factors such as the colour correction. According to the IRAC 
data handbook, the colour correction factor for blackbody spec- 
tra with temperatures between 2 000 and 5 000 K is around 1%, 
except for IRAC 4, where it can reach as much as 2.7% for the 
hotter stars. Since the colour correction is small and we do not 
know a priori the nature and spectrum of our objects, we did not 
apply these corrections. 

2.2.4. Step 4: Creating the point source catalogue 

We finally decided to adopt only the PRF fluxes for the final 
catalogue on the basis of comparison with other catalogues and 
with isochrones in CMDs. With this choice, we find, on the one 
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Fig. 2. Comparison of the fluxes resulting from PRF fitting 
and aperture photometry for the IRAC 1 observations of the 
NGC 6522 field. Top: PRF fitting fluxes and aperture photom- 
etry for the aperture with a radius of 2 pixels; middle: the same 
as top panel, but for the aperture with a radius of 3 pixels; bot- 
tom: relative difference between the measured fluxes (3 pixels 
aperture radius) versus the PRF flux for all sources. The dotted 
lines in the bottom panel mark relative differences of +7%. In the 
upper and middle panels, the dashed line is a linear least squares 
fit to the sources with^ 2 < 5, and the solid horizontal and verti- 
cal lines mark the lower and saturation limits of Table|2]used for 
quality grade labelling. 



hand, slightly better agreement with other catalogues, and on 
the other, somewhat reduced scatter at the faint end of CMDs. 
See Sect. [3] and [4] for a comparison with other catalogues and 
isochrones. 

The last step is the cross-identification, for each field, of the 
extraction tables of all five bands among each other and with 
other reference catalogues. As reference catalogues we adopt the 
Deep Near Infrared Survey (DENIS ; JEpchtein et aL]|1997j i, the 
2 Micron All Sky Survey (2MASS; |Skrutskie et al.||2006} , the 
Midcourse Space Experiment (MSX) catalogue, and the Infra- 
Red Astronomical Satellite (IRAS) point source catalogue. We 
applied the following criteria to include a source in the final cat- 
alogue. We identified individual sources by their RA and Dec 
position and allowed for an error margin of 1'.'6, i.e. slightly 
more than one pixel size, for the cross-identification among our 
five bands. Sources that are detected in at least two of our five 
bands are included in any case in the final catalogue. A cross- 
identification is made for these sources with a search radius of 
3'.'0 with DENIS and 2MASS, and 3070 for identification with 
an MSX or IRAS counterpart. The information about these coun- 
terparts is included in the catalogue. Sources that are detected in 
only one of our five bands are only included in the final cata- 
logue if they also have a counterpart in at least one of the ref- 
erence catalogues, with the same search radii applied as for the 
sources with two or more Spitzer detections. This procedure was 
followed to include as many of the sources in the final catalogue 
that were only detected in MIPS 24. These might very well be 
real objects, albeit extremely red (i.e. not detected with IRAC). 
Additionally, by following this procedure, more of the sources 
that are located in the areas that are not sampled by all five bands 
will be included in the catalogue. These areas are also included 
in the final catalogue. The least reliable sources in our catalogue 
are those with detection in only one IRAC band and no MIPS 
detection, and with (i) either no 2MASS and DENIS counterpart 
(depending on the field, between 3.1% and 6.7% of the sources), 
or (ii) a DENIS and 2MASS counterpart at a distance between 
1 .6 and 3 arcseconds (depending on the field, between 0.4% and 
1.0% of the sources). 

Table [5] gives ten lines of the band-merged catalogue of the 
Bulge N 1 field as an example. The entire table will be made 
available for download from CDS. A portion is shown here for 
guidance regarding its form and content. 

The columns of our point source catalogue are explained as 
follows. 

Column 1: Source identification, IAU-conform identifier. 
Column 2: Right ascension in degrees (J2000). 
Column 3: Declination in degrees (J2000). 
Column 4: PRF flux in IRAC 1 (in yuJy), set to "-9.99999e+99" 
if not detected. 

Column 5: PRF flux in IRAC 2 (in //Jy), set to "-9.99999e+99" 
if not detected. 

Column 6: PRF flux in IRAC 3 (in yuJy), set to "-9.99999e+99" 
if not detected. 

Column 7: PRF flux in IRAC 4 (in yuJy), set to "-9.99999e+99" 
if not detected. 

Column 8: PRF flux in MIPS 24 (in fiiy), set to 

"-9.99999e+99" if not detected. 
Column 9: Quality flags of the IRAC1, 2, 3, 4, and MIPS 24 

channels (see Section [2. 2. 2\ , set to "— " if not detected. 
Column 10: Uncertainty of the PRF flux in IRAC 1 (in /Jy), set 

to "-9.99999e+99" if not detected. 
Column 1 1 : Uncertainty of the PRF flux in IRAC 2 (in //Jy), set 

to "-9.99999e+99" if not detected. 
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Column 12: Uncertainty of the PRF flux in IRAC 3 (in yuJy), set 
to "-9.99999e+99" if not detected. 

Column 13: Uncertainty of the PRF flux in IRAC 4 (in ,uJy), set 
to "-9.99999e+99" if not detected. 

Column 14: Uncertainty of the PRF flux in MIPS 24 (in yuJy), 
set to "-9.99999e+99" if not detected. 

Column 15: Observation flag. The first four bits state if the posi- 
tion of the source was in the field of view of the IRAC bands 
("1") or not ("0"), the fifth bit is for MIPS 24. 

Column 16: Distance to closest DENIS source (arcsec; only 

< 3'.'0), set to "-" if no counterpart. 

Column 17: DENIS /-band magnitude, set to "-" if no counter- 
part or the counterpart has no /-band measurement. 

Column 18: DENIS /-band magnitude, set to "-" if no counter- 
part or the counterpart has no /-band measurement. 

Column 19: DENIS ZT-band magnitude, set to "-" if no coun- 
terpart or the counterpart has no ZT-band measurement. 

Column 20: Distance to closest 2MASS source (arcsec; only 

< 3'.'0), set to "-" if no counterpart. Only sources with qual- 
ity flags A, B, C, or D in at least one 2MASS filter have been 
considered. 

Column 21: 2MASS /-band magnitude, set to "-" if no coun- 
terpart within 3'.'0 was found, or "-9.999" if the counter- 
part has no measured /-band magnitude available or too low 
quality flag. 

Column 22: 2MASS //-band magnitude, set to "-" if no coun- 
terpart within 3'.'0 was found, or "-9.999" if the counter- 
part has no measured //-band magnitude available or too low 
quality flag. 

Column 23: 2MASS ZT-band magnitude, set to "-" if no coun- 
terpart within 3'.'0 was found, or "-9.999" if the counter- 
part has no measured ZT-band magnitude available or too low 
quality flag. 

Column 24: 2MASS quality flags, set to "-" if no counterpart 
or too low quality flag. 

Column 25: Distance to closest IRAS source (arcsec; only 

< 30"), set to "-" if no counterpart. Only sources with qual- 
ity flags 2 or 3 in at least one IRAS band have been consid- 
ered. 

Column 26: IRAS 12yum flux (Jy), set to "-" if no counterpart, 
and to "-9.9e+99" if the counterpart has no measured flux 
in the 12/im band or too low quality flag. 

Column 27: IRAS 12/im quality flag, set to "-" if no counter- 
part or too low quality flag. 

Column 28: IRAS 25 yum flux (Jy), set to "-" if no counterpart, 
and to "-9.9e+99" if the counterpart has no measured flux 
in the 25 fim band or too low quality flag. 

Column 29: IRAS 25 fim quality flag, set to "-" if no counter- 
part or too low quality flag. 

Column 30: IRAS 60yum flux (Jy), set to "-" if no counterpart, 
and to "-9.9e+99" if the counterpart has no measured flux 
in the 60 fim band or too low quality flag. 

Column 3 1 : IRAS 60 fim quality flag, set to "-" if no counter- 
part or too low quality flag. 

Column 32: IRAS 100 fim flux (Jy), set to "-" if no counterpart, 
and to "-9.9e+99" if the counterpart has no measured flux 
in the 100 /im band or too low quality flag. 

Column 33: IRAS 100 //m quality flag, set to "-" if no counter- 
part or too low quality flag. 



Column 34: Distance to closest MSX source (arcsec; only 
< 30"), set to "-" if no counterpart. Only sources with qual- 
ity flags 2, 3, or 4 in at least one MSX band have been con- 
sidered. 

Column 35: MSX Bl band flux (Jy), set to "-" if no counter- 
part, and to "-9.999e+99" if the counterpart has no mea- 
sured flux in the Bl band or too low Bl band quality flag. 

Column 36: MSX Bl band quality flag, set to "-" if no coun- 
terpart or too low quality flag. 

Column 37: MSX B2 band flux (Jy), set to "-" if no counter- 
part, and to "-9.999e+99" if the counterpart has no mea- 
sured flux in the B2 band or too low B2 band quality flag. 

Column 38: MSX B2 band quality flag, set to "-" if no coun- 
terpart or too low quality flag. 

Column 39: MSX A band flux (Jy), set to "-" if no counterpart, 
and to "-9.999e+99" if the counterpart has no measured flux 
in the A band or too low A band quality flag. 

Column 40: MSX A band quality flag, set to "-" if no counter- 
part or too low quality flag. 

Column 41: MSX C band flux (Jy), set to "-" if no counterpart, 
and to "-9.999e+99" if the counterpart has no measured flux 
in the C band or too low C band quality flag. 

Column 42: MSX C band quality flag, set to "-" if no counter- 
part or too low quality flag. 

Column 43: MSX D band flux (Jy), set to "-" if no counterpart, 
and to "-9.999e+99" if the counterpart has no measured flux 
in the D band or too low D band quality flag. 

Column 44: MSX D band quality flag, set to "-" if no counter- 
part or too low quality flag. 

Column 45: MSX E band flux (Jy), set to "-" if no counterpart, 
and to "-9.999e+99" if the counterpart has no measured flux 
in the E band or too low E band quality flag. 

Column 46: MSX E band quality flag, set to "-" if no counter- 
part or too low quality flag. 
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3. Comparison with other missions and catalogues 

Our observations overlap with a number of other mid-IR sur- 
veys. These are 



L. 



ISOGAL: Survey with the ISO Camera (ISOCAM) on-board 
the ISO satellite, combined with DENIS IJK S photometry. 
Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
II (GLIMPSE-II), a Spitzer/IRAC survey of the area +1° 
around the Galactic plane and the Galactic centre (Benjamin 
|et al.|2003) . 

GALactic CENtre (GALCEN), a Spitzer/IRAC survey of the 
1°4 X 2°0 of the Galaxy ( |Ramfrez et al.|2008) . 



inner 



4. A Spitzer/MlPS survey of the inner - 
( |Hinz et al.|2008) . 

3.1. Comparison with ISOGAL 



1.5x8.0 of the Galaxy 



ISOGAL ( |Omont et aLl|20"03) is a survey with the ISOCAM 
instrument on-board the ISO satellite in two bands at 7 and 
15 pm, combined with IJKs photometry from the DENIS project 
( |Epchtein et aT1[l997] l. The 7 pm LW2 band of ISOCAM (~ 
5.0-8.5 pm, central wavelength 6.75 pm) overlaps in wavelength 
with the IRAC4 channel (~ 6.4 - 9.4 pm, nominal wavelength 
7.844jum; |Hora et al.|2008| >. For IRAC4, we have overlap with 
ISOGAL in the fields Bulge 2, Bulge 4, N 1, and NGC 6522. The 
number of sources in common (with less than 270 positional 
offset) is 163 for Bulge 2, 335 for Bulge 4, 288 for Nl, and 
264 in NGC 6522. We present here a brief comparison for the 
NGC 6522 field, because the ISO 7 pm photometry goes deepest 
in that field (least crowding and diffuse background emission). 
The results for the other fields are very similar. 

Figure[3]shows the magnitude difference between the Spitzer 
IRAC4 band and the ISO LW2 band for the NGC 6522 field. 
The data points do not scatter randomly around zero. Rather, 
faint sources tend to be brighter in the ISO 7 pm band, whereas 
bright sources tend to be brighter in the IRAC 4 band. The rela- 
tion found to convert between the two bands is 



magisoi = 0.866 x mag IRA c4 + 1.024. 



(2) 



The slopes and zero points of this linear fit are similar for all 
fields, with a cross-over {Spitzer IRAC 4 equal to ISO 7 pm mag- 
nitude) between 6™0 and 7™7. The trend might be a reflection of 
a dependence of the strength of the SiO fundamental band at 
~ 7.5 pm and/or a water band at ~ 5.5 pm on the brightness of 
the star. At least the SiO first overtone band at ~ 4.0 pm has 
been found to decrease in strength from semi-regular variables 
to Mira-like variables (Aringer et al. 1999). On top of that, in- 
strumental effects in one or both missions cannot be excluded. 



3.2. Comparison with GLIMPSE-II 

The Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 

2003 [ l, is a Spitzer Legacy 



(GLIMPSE-II; Benjamin et al. 



Science Programme of most of the inner Galactic disk using the 
IRAC instrument. Our programme has overlaps the GLIMPSE- 
II survey in the fields Bulge 1, Bulge 2, Bulge 4, and Bulge N 1. 
A comparison with this large data set is especially interesting be- 
cause the same instrument and the same filter bands have been 
used, though the observing strategies and data reduction tech- 
niques are different. 

In Fig. |4] we show the magnitude differences between our 
and the GLIMPSE-II catalogue for all four IRAC channels for 
the Bulge N 1 field, since there is a large number of sources 




Fig. 3. Magnitude difference between the Spitzer IRAC 4 band 
and the ISO 7 pm LW2 band for 264 sources in common in the 
NGC 6522 field. The dotted line marks the identity relation, the 
dashed line is a linear least squares fit through the data points, 
excluding those sources fainter than 9™3 in either of the data sets 
to avoid Malmquist bias. 



in common in this field (31 191). These plots look very similar 
for the other fields with overlap. In general, the agreement be- 
tween our and the GLIMPSE-II catalogue is good in the bright 
to medium brightness range. However, at the faint end, start- 
ing at around 10™ - 11™, we notice a strong trend such that the 
magnitudes become increasingly fainter in our catalogue than in 
the GLIMPSE-II catalogue. This is not only true for the flag B 
sources (most of the red dots at faint magnitudes in Fig.|4|, but 
also for flag A sources fainter than this magnitude limit, and for 
all filters. The same trend is revealed in a comp arisorp] be tween 
the GLIMPSE-II catalogue and the one of |Rarmrez et al.| ( |2008] 
see next section). According to B. Babler of the GLIMPSE-II 
team (private communication, 2008), this trend at the faint mag- 
nitude limit arises because of GLIMPSE-II's limitation to "sin- 
gle frame" photometry. GLIMPSE-II does photometry only on 
single bed frames instead of mosaics, because each patch of the 
sky is observed only twice. This strategy may cause a Malmquist 
bias, such that GLIMPSE-II magnitudes will be increasingly too 
bright the closer a source is to the faint limit. Thus, for sources 
at the faint limit, our magnitudes are likely more accurate than 
GLIMPSE-II magnitudes. Finally, we also notice some satura- 
tion effects above ~ 6™ 5 in IRAC 1, and a few sources that are 
significantly fainter in the GLIMPSE-II catalogue than in our 
catalogue that can be found at all magnitudes. 

Due to the large number of sources in common with the 
GLIMPSE-II survey, we are able to check whether or not the 
combined errors are of realistic magnitude. To do so, we inspect 
the sigma-factor, which is also discussed in the aforementioned 
comparison between the GLIMPSE-II and Ramirez et al. (2008) 
catalogues: 



our mag - GLIMPSE mag 
V(our error) 2 + (GLIMPSE error) 2 



(3) 



If the combined errors are close to the magnitude differences, the 
sigma-factor will have a Gaussian distribution of width 1 . If the 
errors are underestimated, the distribution will be broader than 
that and narrower if the error estimates are too large. Figure [5] 
displays the sigma distribution for all four IRAC channels for the 
field Bulge N 1, along with a Gaussian fit to the data. Also plot- 
ted is a Gaussian for a "perfect" distribution with the same area 
under the graph, i.e. when magnitude differences and combined 
error are on average of the same size (cr - 1). The source selec- 
tion was restricted to the magnitude range 7™0- 10™0 to avoid the 
tail of faint sources with a Malmquist bias (see above). We see 



http://www.astro.wisc.edu/sirtf/glm2_galcen_comparison.pdf 
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Fig. 4. Magnitude differences between our and the GLIMPSE-II catalogue for the sources in common in the N 1 field. Sources with 
quality flag A in our catalogue are represented by black dots, all other quality flags by red dots. 



that the distributions are centred on zero, but there is a surplus 
of negative sigma values, because of the sources with negative 
magnitude difference, as noted above. The actual distributions 
for IRAC 1 and 2 are only slightly broader than what would 
be expected if the errors were correctly estimated. This means 
that the combined errors for these channels are only slightly un- 
derestimated. For IRAC 3 and 4, the distributions are definitely 
broader than in the ideal case. For these channels we have to 
assume that the combined errors are underestimated by factor of 
1 .8 to 2.2. Indeed, our error estimates are on average smaller than 
the ones estimated by GLIMPSE-II, but this is mostly due to the 
longer exposure time per pixel in our survey, and we are doing 
the photometry on the mosaicked images. If we used two times 
the error as given by GLIMPSE-II in Eq. [3] the sigma distribu- 
tion would have the perfect width of 1.0 for IRAC 1, would be 
slightly too narrow for IRAC 2 (too large errors), and would still 
be too broad for IRAC 3 and 4 (too small errors). We therefore 
conclude that our error estimates for IRAC 1 and 2 are proba- 
bly only slightly too small, and too small by a somewhat larger 
amount for IRAC 3 and 4. The GLIMPSE-II errors, on the other 
hand, are probably correctly estimated for IRAC 1, slightly too 
large for IRAC 2, and slightly underestimated for IRAC 3 and 
4. Table [6] summarises the width, centre, and peak values of the 
Gaussian fits to the sigma distributions of Fig. [5] 



Table 6. Width, centre, and peak values of the sigma distribu- 
tions. 



Band 


o~ 


Centre 


Peak 


IRAC1 


1.42 


-0.22 


241.5 


IRAC 2 


1.27 


0.17 


159.0 


IRAC 3 


1.80 


-0.24 


202.7 


IRAC 4 


2.15 


-0.17 


98.6 



3.3. Comparison with Ramirez et al. '2008) 



The overlap with the GALCEN survey of |Ranurez et al.| (2008) 
is limited to the Bulge N 1 field. The reason for this is that 
the Spitzer time allocation committee tried to avoid redundant 
observations of the same fields of sky within GO programmes 
(GLIMPSE-II, on the other hand, is a legacy survey). A compar- 
ison between our own and the GALCEN observations is useful, 
because for our N 2 field we will have to use the IRAC observa- 
tions of GALCEN. Nevertheless, because IRAC channels 1 and 
3 are seeing a patch of the sky neighbouring that of channels 2 
and 4, we got a small overlap with Ramirez et al.| ( [20 08) in our 
Bulge N 1 field in order to fully sample that field in all bands. 
Adopting a stringent 170 search radius, there are 5172 sources 
in common in any one of the four IRAC bands. In the Bulge N 2 
field, 337 sources detected with MIPS 24 have a counterpart in 
Rami rez et al.| ( |2008) , adopting a less stringent 370 search radius 
because of the less precise MIPS coordinates. 
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Fig. 5. Distribution of sigma-factors (Eq.^ for all IRAC channels in the field N 1 in the magnitude range 7™0 - 10?0. The solid line 
is a Gaussian fit to the distribution. The dotted line indicates a Gaussian of width 1 and the same area under the graph, the ideal case 
of correctly estimated errors. 



Figure [6] shows the magnitude differences found for the 
sources in common with Ramirez et al. (2008 ) in the Bulge N 1 
field. The agreement in the IRAC 1 and 2 bands is very good. In 
the IRAC 3 band, the GALCEN magnitudes (of flag A sources in 
the medium brightness range 7^0 - lOTO) are on average ~ O^OS 
fainter than our magnitudes. In the IRAC 4 band, a "knee" ap- 
pears for sources brighter than 8™2. Since the same effect is 
found in the aforementioned comparison between GLIMPSE- 
II and GALCEN data, the origin of this "knee" most probably 
comes from an underestimation of fluxes for sources brighter 
than ~ 8™2 in IRAC 4 by GALCEN. We will have to con- 
sider these differences between our IRAC photometry and that 
of Ramrre z~et al.| ([2008 ) when comparing the results for Bulge 
N 2 to those of other fields. 

The differences between our reduction and that of GALCEN, 
on the one hand, and GLIMPSE-II, on the other, should be noted. 
While in our reduction the background is determined for each 
pixel in a 45 x 45 pixel box centred on that pixel and subtracted 
prior to the photometric measurements, GALCEN uses the local 
background determined by PRF fitting and subtracts it from the 
corresponding aperture flux of each detected source. For the fit- 
ting, PRFs pre-defined by the SSC are used. GALCEN uses the 
aperture corrections as recommended by the SSC, but adopts the 
PRF fluxes for their final catalogue (except for sources where 
the ratio between PRF flux and aperture flux exceeds 1.5; for 
those sources the aperture flux is adopted). We adopted the PRF 



flux in our catalogue, derived from the most recent mean PRFs 
provided by the SSC. GLIMPSE-II, on the other hand, performs 
photometry on individual bed frames, not on mosaics. A com- 
bination of PRF fitting and aperture photometry called "tweak- 
ing']^] is used to measure the flux of a source. PRFs dynam- 
ically determined from each individual frame are fitted to the 
sources and subtracted from the image. On the residual image, 
aperture photometry was performed around the positions where 
PRFs had been subtracted in the previous step. If the resultant 
aperture flux is substantially positive or negative, the source has 
been under- or over-subtracted. This residual aperture flux is 
then subtracted from or added to the PRF flux to compensate for 
the under- or over-subtraction. However, tweaking was applied 
rather sparingly and only in the IRAC 1 and IRAC 2 bands of the 
GLIMPSE-II observations. These differences in reduction have 
to be kept in mind, just as for the comparison between GALCEN 
and GLIMPSE-II. 



3.4. Comparison of MIPS 24 measurements with Hinz et al. 



(2008) 

A catalogue of Spitzer/MJPS 24 fim sources towards the region 
of ~ l°.5 X 8°0 around the Galactic centre has recently been 
presented by Hinz et al. (2008 1. With this survey, we have 371 



2 http://www.astro.wisc.edU/sirtf/glimpse_photometry_vl.0.pdf 
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Fig. 6. Magnitude differences between our photometry and that of 
as in Fig. [4] 



Ramirez et al. 



( 2008 i for the small overlap in field N 1 . Symbols 



sources in common (i.e. unique counterparts within a one phys- 
ical pixel = 2'.'55 search radius) in our Bulge 4 field, 280 in the 
Bulge N 1 field, and 129 in the Bulge N2 field. Thanks to the 
longer exposure time of our observations (33 1 s per pixel com- 
pared to 15 s per pixel in the fast scan mode used by Hinz et al.), 
our catalogue reaches much fainter flux levels. 

For the N 1 field, we only have a small overlap with the sur- 
vey of Hinz et al. ( 2008) . In the N2 field there is strong diffuse 
galactic emission, and a small fraction of the sources might actu- 
ally be false detections. We thus concentrate on the Bulge 4 field 
for a comparison with the |Hinz et al.| ( [2008] ) catalogue. This is 
also the field with the most sources in common. In Fig. [7] we 
show the magnitude differences between our measurements and 
that of Hinz et al. ( 2008) as a function of the source magnitude in 
our catalogue. There might be a general offset of ~ 07*15 in the 
sense that our magnitudes are on average fainter, but the scatter 
is considerable. At least there seems to be no trend at the faint 
magnitude end. 

As for the comparison with GLIMPSE-II, we also checked 
the error estimates of the MIPS 24 data by investigating the dis- 
tribution of sigma-factors (Eq. |3j. The result of this exercise is 
that the combined errors are too small by a factor of six or more. 
Even adopting twice the error estimate of Hinz et aT] (2008), 
which is somewhat larger than ours, instead of a combination 
of the errors, gives a sigma distribution that is wider by a fac- 
tor of four than what is expected from a correct estimate. We 
thus conclude that the errors are still largely underestimated in 



both catalogues, and that the uncertainty due to crowding and 
strong background radiations renders a precise flux determina- 
tion in the vicinity of the Galactic centre impossible. Some part 
of the found differences, however, may be explained by a real 
variability, since many of the bright sources are expected to be 
AGB variables in the bulge (Glass et a l.|2009) . At a wavelength 
of 20 /mi, a full amplitude of up to 07*7 is found for C-rich pul- 
sating AGB variables in the solar neighbourhood (see e.g. Fig. 7 
of |Le Bertre|1992l >. 

4. Colour-magnitude diagrams and comparison 
with theoretical isochrones 

A final check on the quality of our data reduction are colour- 
magnitude diagrams (CMDs) in combination with theoretical 
isochrones. If the flux measurement has no systematic error (and 
also the isochrones are producing realistic colours), the observed 
giant branch will be described well by the isochrone. This check 
is restricted to the IRAC bands, since for many sources the 
MIPS 24 band will be affected by dust emission, which is very 
difficult to model accurately in isochrones. We chose to use the 
most recent isochrones from Marigo et al. ( 2008[ l. A web-form 
can be used for computing these isochrones with different pa- 
rameters^] We computed an isochrone with an age of lOGyrs, 
metallicity Z = 0.019 (solar metallicity), and no dust formation. 



http://stev.oapd.inaf.it/cgi-bin/cmd_2.! 
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Fig. 7. Magnitude differences between our MIPS 24 magnitudes 
and |HinzetaT1 ( |2008l l, for our Bulge 4 field. 



Table 7. Coverage of our fields by different extinction maps of 
the Galactic centre region. 



Field 


Ref. 


median A v 


Bulge 1 


(D,3 


4.38 


Bulge 2 


2,3 


2.45 


Bulge 3 


(1),3 


2.00 


Bulge 4 


1, 3 


15.75 


Bulge N 1 


1,3 


7.25 


Bulge N 2 


1,3 


20.06 


NGC 6522 


1, 2, 3 


1.13 



Refere nces: l: |Schultheis"eTaTJ ( [T999t ; 2:[Sumi]([2004j; 3: |Dutra et al] 
(2003 1. Numbers in brackets mean that the field is not entirely covered 
by the respective map, but to a significant fraction. 



In the mid-IR, the isochrones are not very sensitive to the pre- 
cise choice of the parameters such as age, metallicity, and dust 
formation. 

To disentangle possible systematic errors in the IRAC chan- 
nels, we want to include a flux measurement that is independent 
of our reduction method in the comparison with the isochrone. 
Since all our fields are covered by 2MASS, which can be re- 
garded as a reliable source of Ks magnitudes, we decided to 
construct dereddened K sfi vs. (K s - IRAC) CMDs. 

Most of the fields towards the Galactic bulge suffer from 
strong extinction. Thus, it is necessary to correct the photometric 
measurements for this extinction before comparing them to the 
theoretical isochrone. Several extinction maps for the Galactic 

(1999) de- 



centre region have been published. Schulfheis et al 
rived the extinction from the shift of the giant branch in the 
vs. (J - Ks,) CMD based on DENIS data. The extinction values 
published in that work are given in the V band. Another map 
of Ay values towards fields in the Galactic bulge is presented 
by |Sumi ( 2004| l. These are based on the extinction measured on 



red clump giants in fields of the OGLE-II survey ( |Udalski et al 
20001. F inally, the map ofJDutra eFaLl ( [2003] l uses the same tech- 
nique as [Schultheis et al. ( 1999) and applies it to 2MASS data. 
That map gives the reddening value in the K$ band. Its spatial 
resolution is lower than that of the other maps, but it is the only 
one that fully covers all of our Spitzer fields. 

Table|7]summarises the coverage of our seven fields by these 
extinction maps, and the median extinction in the V band is 
given. The NGC 6522 field is the only one that is covered by 



all three extinction maps. Since this field also has the lowest ex- 
tinction of all our Spitzer fields so that uncertainties in the ex- 
tinction determination should play a minor role, it serves here as 
a "fiducial" field for the comparison with theoretical isochrones 
in CMDs. We used the average values for the diffuse inter-stellar 
medium determined by Indebetouw et al. ( 2005 i for the extinc- 
tion in the IRAC bands: A K s : A3 6 : A45 : A58 : Ago = 1 : 
0.56 : 0.43 : 0.43 : 0.43, as well as A v : A K , S = 1 : 0.089. 
Using the somewhat different values found by other studies (e.g. 
Flaherty et al. 2007 1 Nishiyama et al. 2009) has only a negligible 
impact on the results. Figure |8J shows CMDs of the NGC 6522 
field involving the 2MASS Ks magnitude and the magnitudes of 
the four IRAC bands. 

The following observations can be made from these CMDs. 
At the faint end, the distribution in colour is quite broad since the 
errors become quite large for these faint sources, and some of the 
objects probably are background galaxies. There are also a small 
number of objects with extremely red or extremely blue colours. 
Possibly, false identifications are involved in some cases (we in- 
clude only sources with a 2MASS counterpart within 2'.'0 search 
radius for the comparison with CMDs). At intermediate bright- 
nesses, the distribution in colour is much narrower, and only 
very few outliers can be found. At the bright end, the number 
of sources far from the isochrone locus increases again. This is 
partly due to saturated sources, e.g. in IRAC 1 (blue "sequence" 
at the bright end), and to dust emission that already shows up 
at IRAC 4 wavelengths. In general we find good agreement be- 
tween the observed and the predicted location of the giant branch 
in the CMDs. In the (K s - IRAC 1) CMD (upper left panel of 
Fig. [8j, the tip of the giant branch in the isochrone is somewhat 
bluer than the observed tip. In the (K s - IRAC 3) CMD (lower 
left panel), the theoretical isochrone is bluer than the observed 
giant branch by about 0?1 over the whole brightness range. In 
this filter, the agreement with the GLIMPSE-II survey is very 
good, but the GALCEN magnitudes are fainter than our magni- 
tudes. Thus, by adopting the GALCEN magnitude scale, the shift 
of the isochrone would be somewhat alleviated, but not all cata- 
logues and the isochrone can be made agree. We do not find any 
trends in magnitude in the IRAC 1 filter in the comparison with 
any of the other surveys, thus we assume that the difference is re- 
lated to a small problem with the isochrone in the (Ks - IRAC l)o 
colour. The same might be suspected for the [Ks - IRAC 3)o 
colour. 

The comparison with the isochrones also allows for a check 
on how deep our Spitzer photometry goes. We find that in all 
IRAC bands, we reach sources down to the beginning of the 
He-core burning (horizontal branch), although probably with a 
reduced detection probability and less accurate photometry. The 
He-core burning phase is covered until its end, hence also the 
whole AGB evolution. However, the sensitivity of the present 
observations is still far too low to reach the beginning of the 
RGB phase, and our investigations will be limited to the brighter 
half of the RGB. The MIPS photometry does not reach down to 
the horizontal arm. The faintest sources that are detected with 
MIPS in the NGC 6522 field are roughly 3™5 fainter than the 
RGB tip of the isochrones. 

As for the colour dependence on the used extinction map, 
we find that for individual sources the de-reddened Ks _ IRAC 
colour may vary by a few 0^0 1 . However, no general trend is ob- 
served when different extinction maps are applied. Thus, though 
the fine structure of the extinction varies somewhat from map 
to map, the average magnitude of the extinction is very similar 
among the maps. 
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Fig. 8. Colour-magnitude diagrams of the field NGC 6522 involving the 2MASS K$ and different IRAC magnitudes. For this version, 
the extinction map of Schultheis et al. ( 1999[ ) has been used. A theoretical isochrone from Marigo et al. (2008 ) is included as a blue 
line (see text for details). Symbols are as in Fig. |4] 



5. Summary and conclusions 

We present a catalogue of Spitzer IRAC/MIPS observations of 
seven fields towards the Galactic bulge, sampling a range of 
galactocentric radii. These observations allow us, amongst other 
things, to study the mass loss of a large and homogeneous sam- 
ple of RGB and AGB stars down to lower luminosities and mass- 
loss rates than previously achieved. 

In this first paper, we present the observations, the data re- 
duction procedure, and comparisons to other mid-IR surveys of 
the Galactic bulge. In general, we find good agreement with 
other surveys. The comparison between the Spitzer IRAC 4 band 
and the ISOGAL LW2 band shows good agreement, but re- 
veals a slight trend with magnitude. GLIMPSE-II magnitudes 
are in good agreement with our magnitudes in the bright-to- 
medium brightness range, but a strong trend is present at the 
faint end. This trend is probably related to a Malmquist bias in 
the GLIMPSE-II data set. The error estimates of GLIMPSE-II 
and our IRAC photometry are reasonable in the IRAC 1 and 2 
bands, but somewhat too small in the IRAC 3 and 4 bands. In the 
comparison with GALCEN, we find on average ~ 0^05 brighter 
magnitudes in the IRAC 3 band. The source of the discrepancy at 
the bright end of the IRAC 4 band is probably not related to our 
catalogue. A comparison with the MIPS 24 catalogue of Hinz 
et al. ( 2008 1 reveals that our magnitudes are probably brighter 
at the level of ~ 0^15, although with larger scatter. We also 



find good agreement between our data and recent isochrones in 
colour-magnitude diagrams for at least three of the four IRAC 
bands. The ~ 0^1 offset from isochrones involving IRAC 3 de- 
serves some more attention. We thus may assume that the ob- 
servations and data reduction are accurate on the level of ~ 0?1 
or better, as well as precise on the level of < 07*1, except for 
faint sources with quality grade B. The science exploitation of 
the data will follow in subsequent papers. 
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Table A.l. Detected and extracted sources. 





detected 


extracted 


SNR<5 


sources outside 


both criteria 


no deblending 


37974 


37471 


380 


174 


51 


pass, deblending 


37974 


37314 


483 


218 


41 


act. deblending 


70474 


43655 


2145 


25885 


1211 


act. & pass, deblend. 


49679 


39003 


1743 


9900 


967 



Appendix A: De-blending or not de-blending? 

We tested the comparability of the fluxes derived by PRF fit- 
ting and those resulting from aperture photometry for the IRAC 
3.6 fim channel and for the field NGC6522 using the two de- 
blending capabilities in APEX. The first method is passive de- 
blending, where already during the detection process multiple 
peaks in a cluster of pixels are found and this information is used 
in the PRF fitting module to simultaneously fit all the peaks. The 
second method is active de-blending, i.e. the number of assumed 
sources is iteratively increased during PRF fitting to decrease the 
final x 1 ■ 

We use the point source probability images in the detection 
process and filter the raw extraction table with the criterion that 
SNR is greater than 5 and that APEX has not flagged the source 
as a point source outside the fitting area. This excludes very 
weak sources, sources covered by only one image, and sources 
whose position is not correctly identified. The numbers of de- 



tected and extracted sources for the runs also shown in Fig. A.l 



are given in Table A. 1 along with the number of sources filtered 
by the two criteria. 

The active de-blending procedure relies on the x 2 values to 
decide whether de-blending is necessary or not, which results 
in heavy de-blending activity, even for sources whose PRF and 
aperture fluxes are consistent. We therefore disabled active de- 
blending in the final data reduction. 

Passive de-blending uses the information from the detection 
process about whether different entries in the detection table be- 
long to the same cluster of pixels above the detection threshold. 
If yes, APEX fits the cluster simultaneously, if not, the single 
sources one after the other. The first effect of switching on pas- 
sive de-blending is that the number of faint sources that lie be- 
low an SNR value of 5, hence the number of sources rejected 
by the first criterion, increases (see Table |AT) i. Since the posi- 
tions of the detected sources are also varied to fit the complete 
cluster, the number of sources whose pointing is outside of the 
tile used for PRF fitting increases, too (Table |A.1[ ). The num- 
ber of extracted sources is thus less than in the first run when 
passive de-blending is not activated. Also brighter sources are 
affected, when there are other stars within 5-10 pixels. Then the 
independent stars belong to the same cluster and are fitted simul- 
taneously, although the PRF of one star does not contribute at all 
to the flux in the small tile of another star used for PRF fitting. 
This again changes the position of the centroid of the central 
bright star in the cluster and also the flux measured by PRF fit- 
ting. Since the displacements are smaller than the pixel size, the 
aperture flux does not change as much as the PRF flux. The result 
is that the bright star moves away from the line of good agree- 



ment in the fpRF-f a p el -t ure -plane (Fig- A.l '> so me scatter around 
this line is increased. The same effect occurs for saturated and 
super-saturated sources with a double peak leading to the ob- 
served spreading of sources in the top-left part of Fig. [A~T| Thus 
the results of the run with passive de-blending are not as reliable 
as the results of the run without any de-blending. 
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Fig. A.l. Comparison of the fluxes resulting from PRF fitting and aperture photometry without de-blending (first row), with passive 
de-blending (second row), with active de-blending (third row), and with both (last row) for the apertures with radii of 2 pixels (left) 
and 3 pixels (right) for the field NGC 6522 and IRAC 1 . 
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Appendix B: Magnitude histograms of the seven 
bulge fields. 
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Fig. B.3. Histogram of magnitudes for the Bulge 3 field. The bin size is (F2. 
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Fig. B.6. Histogram of magnitudes for the NGC 6522 field. The bin size is 0^2. 
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Fig. B.7. Histogram of magnitudes for the Bulge N 2 field. The bin size is CP2. 



